Modeling semiconductor thermal properties. The dispersion role. 
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.Sh ■ Abstract 

We study heat transport in semiconductor nanostructures by solving the Boltzmann Transport 
Equation (BTE) by means of the Discrete Ordinate Method (DOM). Relaxation time and phase 
and group velocitiy spectral dependencies are taken into account. The Holland model of phonon 
, relaxation time is revisited and recalculated from dispersion relations (taken in litterature) in order 

to match bulk silicon and germanium values. This improved model is then used to predict silicon 
CO ■ nanowire and nanofilm thermal properties in both ballistic and mesoscopic regimes. 
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I. INTRODUCTION 



Semiconductors are often used in high technology due to their mechanical (crystal struc- 
ture) and electrical (band gaps) properties. They are found in many domains such as nano- 
junctionsi, transistors^, or in solar panels^, and are even used in medicine^. The increasing 
use of semiconductor micro/nano-structures has brought, nowadays, a good knowledge on 
charge transport in spite of heat transfer effects at very short scales 5 . The scope of this 
paper is to contribute to the charasterisation of semiconductor nanowire and nanofilm ther- 
mal properties. Semiconductor nanowires can be found for example in transistors^ and 
semiconductor nanofilms in solar cells. 

At these low scale, Fourier's law may give an inaccurate description of heat transfer 
through these nanostructures. One has to deal with more fundamental physics^. Sev- 
eral approaches as lattice dynamic simulation^ 1 ^^- or solving the Boltzmann Transport 
Equatior*i2^i^LLIikL£ (BTE) have tried to explain experimental measurements^ 1 ^. In this 
paper, the BTE, which is a statistical physics equation, is solved by means of the discrete 
ordinate method (DOM) which is widely used in thermal radiation2&2i We employ an im- 
proved model for describing phonon conduction by taking into account spectral dependancy, 
found in velocities and relaxation times, using dispersion curves. 

At macroscopic scale, conductive heat transfer is governed by the heat equation: 



where k is the thermal conductivity. This equation can be derived in the framework of 
statistical physics considering that individual heat carriers (here phonons) are submitted to 
a random walk—. In this diffusive regime, where collisions are numerous, the heat flux is 
related to the temperature gradient by Fourier's law : 



When heat carriers do not undergo collisions, they cross the system like photons in thin 
optical layer. Heat carriers are then said to have a ballistic flight. A radiative-like form 
describes the steady state solution : 



dT(r, t) 
dt 



V[WT(r,t)], 



(1) 



kVT(r,t). 



(2) 




(3) 
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When phonon mean free path is of the same order as the system length, then some of them 
undergo collisions, while others cross the structure without any interaction. The system is 
in the mesoscopic regime. 

Theoretically, as long as wave effects are negligible, the BTE applies to the ballistic, 
mesoscopic and diffusive regimes. A key feature of this approach is to introduce a specific 
intensity for phonon^, similar to photon intensity used in radiative transfer- 4 . Thus, the 
energy flux per apparent surface unit at point r in direction u reads: 

Iu, p {r, u) = hun^r, u)^-, (4) 

with h the Planck constant, n w (r, u) the phonon density at position r going in direction 
u and V g the heat carrier group velocity. The BTE collision term is classicaly approached 
by the single relaxation time approximation and, therefore, in steady state, the transport 
resolution reads: 

1° - I 

u ■ Vlu,, p = — , (5) 

where 7° is the intensity at local equilibrium and r w p the relaxation time. 



II. NUMERICAL RESOLUTION 

A. Nanostructure geometry modelling 

We here consider nanofilms and nanowires made of semiconductor material. The wire 
geometry refers to a cylinder whose length (L) is much larger than its diameter (D), L » D. 
Cylindrical coordinates are obviously well adapted to this case. Films may also be viewed 
as cylinders in the limiting case where their thickness L is much lower than the diameter 
D. Therefore, BTE in cylindrical coordinates still applies, it reads, for an axisymmetric 
problem^: 

r dr +4 dz r <90 +W^~W^J> W 
where k WiP = — is the extinction coefficient and where /i, rj and £ (EqnJTJ) are the direction 
cosines of the propagation direction u (fig. [1]). They are linked to polar angles and i/j 
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characterising u in a local frame by: 

ji = cos(0) sin(V'), 

T) = sin(0) sin^), (7) 
£ = cos(^). 

The term ^ represents the angular redistribution and the variation over 9 (figH]) is not 
written since = in an axisymmetric problem. 

Figure [1] shows a typical geometry under consideration. The two ending sections are 
always considered as black surfaces set at different temperatures: T hot on the left side and 
T co id on the right side (Thot > T co id) . They emit phonons with a blackbody intensity at their 
wall temperature (EqnJH]). 

where V p is the phase velocity and Twaii is either Thot or T co u ■ 

The lateral (cylindrical) surface is treated as adiabatic, which amounts to consider pure 
reflection at the wall. In the wire configuration, wall reflection is in general partly specular 
and partly diffusive. In the film case, it is always specular so that there is no lateral boundary 
influence. Condition on the cylinder axis is due to symmetry specular reflection. 



B. Discrete Ordinate Method - DOM 



The DOM is based on the selection of a finite set of propagation directions (s m , 
m — 1, . . . , M) and corresponding weights (w m )^^. Sn quadrature (fig. EJ) is the most 
commonly used. It is constructed with a maximum of symmetry requirements (to avoid 
directional bias) and the weight selection rules tend to preserve the exact values of some 
direction cosine key moments. 

Boundary condition for any discrete direction leaving adiabatic surfaces are expressed as: 

I m {xp) = - V] w m >I m i(x P )\s' m ■ n| 

m' if sj n n<0 

+ (l-p)/ A (x P ), (9) 

where, due to Sn set symmetries, if s m belongs to the quadrature, so does the specularly 
reflected direction s„. p is the ratio of diffuse to specular reflection (p = yields pure 
specular reflection). 
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In this equation, as in the rest of the text, spectral indices u and p are omitted for sake 
of clarity. 

In this work, S$ quadrature is employed. It involves 80 directions, but thank to the 
symmetry only 40 directions are used, with their associate weight doubled. 



C. Numerical procedures 

The DOM is widely used to model radiative heat transfer and has already been employed 
to resolve steady or unsteady BTE for phonons2i^. Following Lathrop's guidelines^, a 
variable weight scheme (Eqn. [10]) is used to relate cell-face intensities to the central value 
l m ,p in a node P, for a given propagation direction m, at a given frequency u and polarization 
p. Considering phonon propagation directions s m represented in figure [31 the discretized 
intensities will be written as: 

j j Im,P Ifn,W 
lm,E — Im.W H 5 

a 

+ m - r h m '" , (10) 

Im,P ~ I m -l,P 



2 ' 



The four principal cardinal directions (N, E, S, and W) refer to the cell nodes surrounding 
P (fig. Indexes m± \ point in the two directions around m in a same latitude (£ = C te ). 

Direction are swapped with increasing values of £ and, at a given latitude, with increasing 
values of \i. In each direction a finite volume integration of the BTE yields: 

^wlm,w + A m _i/ m _i P + Xsl m ,s + X I°(T P ) 

*m,P — " " a 71 > (11) 

with: 

X w = Az (\p, m \— - /i m Arj , 
ArAz 



2 WmC 



VmW m (c -I) - a 



A 5 = |£ m |Ary, 

A c = rpK u ArAz, (12) 
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and 



a.! = 0. 

2 

Interpolation weights a, b and c have the following values, which guarantees positivity of the 
solution while keeping, as much as possible, a 2 nd order accurancy: 



a = max 0.5, 1 



b = max 0.5, 1 



\/i m \r w Az 



\ + 2\U\rpAr - 2 AzAr 
\Cm\r P Ar 



a i i I ' 

m +2 



A + 2\/j, m \r E Az - 2AzAr 

c = max I 0.5, 1 + 



a i i I ' 

m +2 



A + 2|£ m |r P Ar + 2|/i m |r£;A2 



(13) 

As explained by Lemonnier— , at each new latitude, for a first direction m, the value of 
I m -^P, corresponding to r\ — 0, has to be initialised. This particular direction is in (r, z) 
plane. Consequently, the angular redistribution contribution (— is null. Therefore, 
by setting a m _i = and ji m = — a/1 — I m _i P is computed with a and b becoming: 



a = max 0.5, 1 — 



b' = max 0.5, 1 — 



Ao + 2\U\r P Ar + AzAr^l - 

\Cm\rpAr 

\ + (r E + r w )Az^/l-e 



(14) 



D. Energy conservation 

At this step, we can determine spectral intensities at any nanostructure points, for a given 
direction m, at a frequency u and polarisation p. The integration over the solid angle dfl 
and over the spectrum (du and p) yields a total intensity in a point P. To be able to obtain 
Im,Pi temperature T is required and is deduced from the energy conservation(Eqn{T5]). 

| + V-q = 0, (15) 
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where e stands for the integrated phonon energy density, and q, the phonon total flux 
density: 

r^max r 

q(r) = V / / 4,p(r,fi) • QdSldu. (16) 

In steady state, the divergence of the phonon total flux density must be null. Therefore, 
integrating Eqnj5] over all directions and the entire spectrum yields equilibrium: 

/ / Kj^pdnduj = J2 / ^Z{T)dnduj. (17) 

JO J47T JO J4tt 

This relation is then inverted to get point P temperature from local intesity (I Ut p) knowledge 
in all direction, all frequencies and all relevant branches (Eqn JTTl) . 



E. Dispersion relations 

Many studies have been done for calculating semiconductors thermal properties, using 
linear— i22j2I or linear by section 3 ^ 3 -* 3 ^ spectral dispersion. Fewer litteratures are found with 
truly non-linear phonon dispersion^. In this work, we use a polynomial phonon dispersion 3 ^, 
for silicon (FigHj). For each polarisation branch, the wave vector is discretized into 60 
equally spaced values. For germanium, as suggested by Lacroix— , cubic splines are fitted on 
experimental data 36 (Figj5]). Germanium wave vector is equally shared into 300 bands. Some 
assumptions are introduced, semiconductors are assumed as isotropic, where the studied 
propagation direction is along [100] lattice direction, and only bulk acoustical modes are 
used. Optical mode are neglected (their respective branches are ignored). 

Real phonon dispersion curves for semiconductor are non-linear. Therefore, the group 
velocity V g (ui,p) — is dissociated from the phase velocity V p (u,p) = j^, both being depen- 
dent on wave frequency. As a consequence, the blackbody intensity (EqnJH]) has a complex 
spectral behaviour. The analytical integration over the wave frequency of the equation ( fP7j) 
is not available. To overpass this complex behavior, V p (u,p) is oftenly approximated, in 
litterature, to a constant V g (u,p) (linear by section dispersion curves) or to a constant av- 
eraged velocity, whereas, in this work, we compute the equilibrium equation. As for the 
phase velocity, the group velocity spectral variation is introduced in the extinction coeffi- 
cient = — . This last definition shows that the extinction coefficient K WtP is a spectral 
and branch dependent quantity through V g (p,u) and the relaxation time t UjP . 
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Following Matthiessen's rule, r W) p may be expressed as: 

C = ru\uj,p) + t n \uj, P ) + rr\uj,p). (18) 

In litterature, Holland relaxation times 32 are oftenly applyed. Therefore, as a starting point, 
we use, in this work, the same form (Tab. [I]) to define impurities, Normal and Umklapp 
processes^. Boundary relaxation time is not used since boundary conditions have been 
imposed further up. Note that the relaxation times can also be temperature dependent. 

III. NEW PHONON RELAXATION TIMES IN DISPERSIVE SEMICONDUC- 
TORS 

In this work, the solution of the spectral BTE requires group and phase velocitie determi- 
nation, deduced from the chosen dispersion curve fit, and of the relaxation time approxima- 
tion. This latter is done with Holland formulas. A compatibility problem thus occurs since 
Holland's dispersion curve fits and ours are not similar. Therefore, one can wonder about 
the relaxation time constants A, B^, B?, and Btu to input in the numerical simulation. For 
example, the transition between N and U process for transverse phonon is not the same. 
Holland assumes that the cutoff frequency is u% = B h 1 whereas Han and Klemens^ state 
it occurs at K ™ ax . 

Thus, these differences bring a new frequency variation and also a new spectral domain. 
Consequently, setting this phonon dispersion to obtain thermal properties with Holland's 
method yields to an improper answer— (Fig 16] and FigJTj). Therefore, the relaxation times 
have to be fitted to approach semiconductor conductivity with Holland's relaxation time 
forms but using non-linear dispersion curves. 

A. Silicon 

1. Silicon new relaxation times 

Since, factors F (boundary) and A (impurity) are present in all the different conductivities 
hr, kru, and ki defined by Holland^, we adjust them in a first place using Holland's method 
to predict conductivity (hr is the conductivity due to impurities, boundaries and Normal 
transverse collisions, k TU is the conductivity due to impurities, boundaries, and Umklapp 

8 



transverse collisions, and kz accounts with impurities, boundaries, and longitudinal phonons 
collisions). Boundary and impurity interactions are dominant at low temperature^^i^, 
therefore F and A are fitted in this temperature range. 

Around 10K, close to transverse phonon Normal process highest value, factor is fixed. 
At high temperature, the coefficient Btu is set where Umklapp processes for transverse 
phonons prevails. It is unimportant to determine Bt before Btu, or conversely when one 
process is dominant the other is negligible. On the other hand, factor Bl, corresponding to 
longitudinal contribution is fitted at last, hence ki major portion has transversal conduc- 
tivities non negligible. This last coefficient is found around 100K. 

At very high temperatures, optical phonons are not taken into account. Although con- 
ductivity values correspond to experimental data, the model does not describe entirely the 
submitted physics at these temperatures. 

2. Silicon thick films and silicon thick cylinder 

To simulate a bulk using our new parameters (TabHTl). we settle the cylinder with the 
film condition (p = 0) for a thickness of L = 7.16mm. In figure ([6]), for temperatures 
above 100K, one can notice that the results obtained for the film correspond to Holland's 
predictions, based on measurements^ 1 ^. For lower temperatures (T < 100K), the values 
obtained with BTE resolution are above the reference curve. Even if the dimensions are 
huge, Holland is working on a bulk with an equivalent sample size L = 7.16mm, which 
appears in the boundary reflection contribution (r^ 1 = V S /LF). In our geometry, the 
sample size L corresponds more to our diameter D. Therefore, to simulate a wire, the ratio 
of diffuse to specular reflection p is set equal to one. Thus, having D = L = 7.16mm, 
conductivities obtained for T > 100K, as for the film, match with measurements^ 1 ^ (Fig. 
|6]). For T < 100K, the results are now below the curve. 

A first set of calculations has been performed with a diffuse to specular reflection ratio 
p = (film) and has brought conductivities above or equal to the reference curve. A second 
set has been examined with p — 1 (wire/cylinder, D ~ L) and has given results below 
or identical to the curve. It can be then assumed that there exists a ratio of diffuse to 
specular reflection between and 1 which will give conductivities similar to experimental 
data (the geometric parameter F is not use in BTE resolution since the boundary conditions 
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are already made in section III Aj) . 



B. Germanium new relaxation times 

A difficulty is found when we try to fit germanium relaxation time parameters with the 
equations employed by Holland. The obtained results always underestimate conductivity 
(FigJTJ), with L = D = 2.4mm (Holland's equivalent sample length and relaxation param- 
eters) or with L = D = 3.8mm (Asen-Palmer's equivalent sample length and relaxation 
parameters). Germanium thermal property alters with various models and experimental 
data. Varying the sample length modifies conductivity values at low temperature where 
boundary scattering is important. Furthermore, germanium conductivity seems to be sensi- 
ble to the material doping^^^^ which increases considerably conductivity curve highest 
point. Thus, changing the impurety parameter A alters considerably the conductivity around 
10K. Fitting then the three first parameters F, A, and E>t to obtain correct values for low 
temperature undergo a mismatch around 140K to 300K. 

To avoid this trouble, Singh^ proposes different relaxation time forms based on three 
different sections of the dispersion curves. Here, as for silicon, we want to use Holland's 
relaxation time forms (TabU]) to model germanium thermal properties, where only four 
parameters are needed against 10 for Singh. Therefore, we solve directly the BTE to fit our 
parameters. Inquired values are based on Glassbrenner and Slack experimental data^^ 1 ^. 
The cylinder dimensions have the same lengths than their sample (L = 20mm and D = 
4.4mm). 

A first fit is done for T > 100 K, where conductivity is less sensible to boundary collisions 
and to Normal transversal processes. Thus, Btu, Bl and A parameters are the conductivity 
action switchs over 100K. The dominant band for tranverse normal process conductivity 
contribution is under 100K, which permits to change E>t without disturbing our new values 
obtained for T > 100K. 

These new germanium relaxation time parameters are given in Table IIHI The resulting 
conductivities have less than 5% of relative difference with experimental data for T > 100K 
(Fig. 0). Just below 100K, a maximum of 7% of relative difference is seen. Under 60K, 
Glassbrenner's conductivities are always between our film and wire values. One can say that 
there exists a single ratio of diffuse to specular reflection, for each temperature, which will 
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give conductivities similar to these experimental data. Therefore we are able to valid our 
germanium parameters. 

Note that the spectral discretization is finer for germanium than for silicon. The first 
Brillouin zones are of the same order, but silicon wave vector is divided into 60 equal bands 
whereas, for germanium, it is splitted into 300 strips. Using 300 bands for silicon changes 
conductivity values by less than 1%. Therefore calculating over 60 wave frequencies saves 
calculation times while keeping accuracy to an acceptable level. On the other hand, using 
only 60 bands over the spectrum can change by up to 10% of the conductivity values. To 
take a maximum of spectral information, we have choosen a large number (300) of bands. 



IV. HEAT TRANSFER IN SILICON NANOSTRUCTURES 



A. Films 

New parameters (Tab JIB have been fitted to provide correct thermal properties in silicon 
and they are now used to study films with BTE resolution. As said earlier, BTE gives 
the possibility to study ballistic, mesoscopic and diffusive problems. Therefore, a thickness 
variation brings us to study nanofilms to thick films (or bulk). Ballistic phenomena is then 
more expected for nanofilms whereas diffusion occurs in thick films. Another way to observe 
ballistic behavior is to work with low temperatures. In these condition, low frequency modes 
are dominant (due to Bose-Einstein distribution) which means that phonon mean free path 
is larger (due to the relaxation time frequency dependence). On the other hand, at high 
temperatures, high frequency modes are dominant, which amounts to small phonon mean 
free path. Three phonon collisions are then more frequent, which drive the process to a 
diffusive scheme. 

Figure [8] shows silicon surface unit conductance (W/m 2 K) versus its film thickness for 
several temperatures. According to the film assumption, the section perpendicular to z is 
supposed to be infinite. Therefore, conductance (W/K) is not available. In a film, the 
surface unit conductance G" reads: 

G" = = ( 19 ) 

-AT Az y ' 

Furthermore, in diffusive regime, when Fourier's law (Eqnj2]) is matching the heat 

transfer, the conductivity, at a given temperature, is constant. In that case we get 
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ln(G") = — ln(Az) + ln(fc). This result correspond to the right part of FigureQ. 

For small lengths and at a given temperature, the surface unit conductance becomes 
constant when thickness becomes thinner. This scheme is due to ballistic phenomena and 
is viewed on curves left side. On the contrary to Fourier's regime, when G" is constant, the 
conductivity concept is questionnable since it depends on the system size (Eqn JT^jl . 

The curved parts of G" (FigJHJ) are the results obtained in mesoscopic regimes. In this 
scheme, phonon mean free path is of the same order than the film thickness. Some phonons 
act as they were in a diffusive system. On the other hand, many phonons have a purely 
ballistic behavior. Consequently, the thermal conductivity concept is again questionnable. 

Another benefit in using BTE resolution is that temperature fields in nanostructures 
can be obtained. In a diffusive regime, the temperature gradient is constant. In ballistic 
regime, it is the temperature which is mostly constant. Figure [9] shows thermal profiles 
along the z axis for a two-micron-thick-film. We can clearly see that Fourier's regime is 
not reached before 300K, which correspond almost to the right mesoscopic limitviewed on 
Figure (IEJ). Ballistic regime is nearly reached at 10K. Between these two temperatures, 
mesoscopic regime prevails. This result confirms that thermal conductivity terms have to 
be used with care. It seems that in silicon films, it is not appropriate to use this concept 
below the micrometer scale. 

In comparison with a simple model^, the surface unit conductance, calculated with a 
single linear group velocity, is overestimated compared to our model (FigllO]). At 300K, 
some points are 30% higher than ours for the film simulation. These results show how 
important is to take into account the spectral dependency. Consequently phonon mean free 
paths (£) are also spectrally dependent since £ = V g {uj,p)T{ijj,p). 

Note that for a thickness beneath 20.10 _9 m, there are less than 40 primitive cells across 
the nanofilm. In that case, bulk dispersion properties are, in principle, no longer valid. 
Correspondant results for these thicknesses are given (FigjS]) only to show the numerical 
behavior of our code despite a physical concordance. 

B. Wires 

Silicon new relaxation time parameters are used to determine nanowire thermal proper- 
ties. In comparison with nanofilms, speaking about thermal conductivity in nanowire could 



12 



be taken out of sense. Nevertheless, representing the surface unit conductance versus the 
length do not present much interest for wires, since the thermal property changes also with 
its diameter. Therefore, we plot the temperature field to see, with the thermal gradient 
along the z axis, in which regime, ballistic, mesoscopic or diffusive, the nanowire, at a given 
temperature, is. For a 2micron long nanowire whose diameter is 115nm, Fourier's law ap- 
plies on a large temperature band (Fig JTTj) . Using thermal conductivity concept to describes 
nanowire thermal properties seems then possible. 

A temperature rise leads to Umklapp process domination which is a resistive process and 
then yields a diffusive regime. Similarly, a diameter reduction increases phonon collisions 
with the border, also a resistive process, and favours a diffusive regime. Therefore, ther- 
mal conductivity as silicon nanowire thermal property is suitable for larger temperature or 
smaller diameter. 

For a given diameter, a shorter wire gives more possibility to have phonon mean free path 
of the same order as the wire length. Low temperature can also accentuate this character 
as phonon traveling length grows up. Other results have been obtained with D = 500nm 
and L = 2/zm, or with D = 37nm and L = 150nm. Their temperature fields have the same 
shape than those shown in figure fflTl) . This tends to consider that thermal conductivity is a 
valid concept in nanowires. However mesoscopic and ballistic regime can always be obtained 
in any nanostructures by reducing hugely the dimensions or/and the temperature. 

We have calculated the conductivity for 2pm length nanowires at diffenrent temperatures 
and diameters (Fig JT2l . We retrieve a same profile than analytical, experimental and 
Monte Carlo 47 data. However, a better agreement is found at high temperature for large 
diameters and at low temperature for small diameters. The boundary has been set here as 
purely isotropically scattering surface. 

As seen in section IIII A 2\ a diffuse to specular reflection ratio p smaller than one is 
expected for low temperatures. A fit on p (Tab llVI) . is then done on our 2/im long and 115nm 
of diameter nanowire. It is used for other diameters (Fig JT3l) . Conductivity matches better 
at low temperature for large diameters, with an adjusted p than with p fixed at 1. Note that 
the ratio decreases with temperature. When temperature goes down, low frequency phonons 
are dominant, and their wavelength average increases (due to Bose-Einstein distribution). 
Consequently, phonons will scatter less with the borders since some of their wavelengths will 
be larger than the nanowire roughness. Therefore, the ratio of specular to diffuse reflection 
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goes down with temperature. 

On the other hand, for small diameters, the conductivity obtained does not match well 
experimental data. This can be explained by the fact that with a small diameter, bulk dis- 
persion data are no longer pertinent. At small scale, phonon modes become discrete whereas 
bulk dispersion is continue. The contribution of non-existing modes probably overestimates 
heat transfer. 

For a 37nm diameter nanowire with a ratio of specular to diffuse reflection set equal to 1, 
the difference between our model predictions and Yang's simple model^ are below 7%, but 
with p = 0.8, ot reaches 13%. For every thin wire, boundary scattering has an important 
resistive role which favours diffusive regime. If we further change the boundary conditions, 
a larger difference is obtained between our model and Yang's. Therefore, neglecting the 
spectral dependency can yield a rough answer. 



V. CONCLUSION 



This paper presents the resolution of Boltzmann Transport Equation with the Discrete 
Ordinate Method in semiconductors. Spectral and temperature dependencies have been 
taken into account to define thermal properties. The polynomial function and the cubic 
spline approximations, representing frequencies u versus wave vector, bring us closer to 
reality. Different relaxation time models, which are made to take into account spectral 
dependency, are proposed in litterature. We have focused on the same relaxation time forms 
applied by Holland, which is oftenly used in literature. Thus, new relaxation times have 
been found to balance non-linear dispersion relations. Sensibilities over parameters defining 
semiconductor relaxation times are different for silicon and germanium. Therefore, for each 
semiconductor, a strategy has been developped to find its own parameters. A comparison 
with a simple model has put in forward the spectral dependency problem. Differences over 
30% can be obtained. It seems that simple models are correctly fitted for specific samples 
or temperatures. Our new parameters, associated to this paper model, permit to describe 
correctly germanium and silicon thermal properties over a large temperature band and over a 
large spatial scale. However, this method does not treat optical phonons and therefore it does 
not represent properly high temperature phenomena. Futhermore, the dispersion relations 
are taken from the bulk modes. Working on very short scales can undergo incorrect answers 
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since phonon dispersion changes. 

Silicon nanofilm surface unit conductance as silicon nanofilm field temperature have 
shown that thermal conductivity is not always a relevant quantity, since diffusive regime 
is only reached, at 1500K, for one micron thick. Even if optical phonons are not described, 
it had been viewed that at lower temperatures diffusive regime is only reached for thick films. 
On the contrary, we have seen that thermal conductivity is a relevant quantity in nanowires 
due to phonon collision with borders. In extremely small nanostructures, we have noticed 
that our treatment can be inproper. Phonon wave behaviour should probably be taken into 
account. Note that the method developed here can easily be generalized for transient heat 
transfer in order to study heat pulse propagation in nanostructures^. 
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FIG. 1: (Color online) Geometry representation 




FIG. 2: (Color online) 56 quadrature (first octant) 



19 




FIG. 3: Four principle cardinal directions surrounding point P with a given propagation direction 
s m . W and S are the known values. Their positions are fixed in opposition to s m propagation 
direction. 




FIG. 4: (Color online) Silicon dispersion relation in direction (100) given by E.Pop 35 
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FIG. 5: (Color online) Germanium dispersion relation in direction (100) given by G.Nilsson 36 




FIG. 6: (Color online) This study is done in a silicon structure, where the dimension parameters 
equal to 7.16mm. The solid line is the conductivity obtain with Holland's method, relaxation 
times and phonon dispersion. The dot line is the conductivity obtain with Holland's method and 
relaxation times but with Pop's phonon dispersion. Dots (•) and (o) are the conductivity obtain 
with BTE resolution but with our new relaxation times and Pop's phonon dispersion. The wire is 
corresponding to p = 1 represented by (•) and the film to p = drawn with (o) (Eqnj9|). Dots (+) 
are Glassbrenner experimental conductivities. 
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FIG. 7: (Color online) This study is done in a Germanium structure. The solid line is the conduc- 
tivity obtained with Holland's method (relaxation times and phonon dispersion), with a sample 
length of L = 2.4mm. The dot line is the conductivity obtained with Holland's method but with 
Asen relaxation time parameters and its dispersion. Asen sample length is L = 3.8mm. The dash 
dot line is the conductivity resulting from Holland's method with Nilsson's phonons dispersion 
{L = 2.4mm). The hyphen line is the conductivity calculated with Asen parameters on Holland 
relaxation time forms with Nilsson's dispersion (L = 3.8mm). Dots (•) and (o) are conductivities 
obtained with BTE resolution but with our new relaxation times and Nilsson's phonon dispersion, 
on a sample of L = 20mm and D = 4.4mm. The wire is corresponding to p = 1 represented with 
(•) and the film to p = drawn with (o) (Eqnj9|). Cross dots (+) are Glassbrenner experimental 
data on a sample of L = 20mm and D = 4.4mm. 
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FIG. 8: (Color online) Film surface unit conductance versus thickness at different temperature. 
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FIG. 9: (Color online) Axis temperature field in a 2micron silicon film. 
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FIG. 10: (Color online) Surface unit conductance comparison between a simple model(Yang) and 
our model(Terris), of films set at 300K. 
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FIG. 11: (Color online) Axis temperature field in a 2fim length and 115nm diameter nanowire. 
The ratio of diffuse to specular reflection is equal to 1. 
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FIG. 12: (Color online) Nanowire themal conductivity. Comparison with experimental results of 
Li et al^, analytical data of Chantrenne et al<^ and Monte Carlo simulations of Lacroix et al4£ 
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FIG. 13: (Color online) Nanowire themal conductivity. Comparison with experimental results of 
Li et al.— , analytical datats of Chantrenne et al.— and Monte Carlo simulations of Lacroix et al.— 



Impureties 


-r 1 


(Rayleigh's model) 


= Au 4 


Longitudinal 


Normal+Umklapp 


polarization 


= B l lu 2 T 3 




Normal 


Transversal 


= BtcoT 4 if < K < K Y X 


polarisation 


Umklapp 

= B u fl^\ ii^<K< K max 



TABLE I: Relaxation time forms. (K is the wave vector, with K max corresponding to the first 
Brillouin's zone, and ks is Boltzmann's constant. K max = where the pure lattice parameter 
a Q = 543pm for silico n 37 ' 38 and a Q = 565pm for germanium 3 -^.) 
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Holland 


Present work 


F 


0.8 


0.66 


A (s 3 ) 


1.32 10" 45 


1.498 10" 45 


B T (K- 3 ) 


9.3 10~ 13 


8.708 10~ 13 


Btu (s) 


5.5 1(T 18 


2.89 10~ 18 


B l (K- 3 ) 


2.0 10~ 24 


1.18 10~ 24 



TABLE II: Parameters for silicon relaxation times 





Holland 


Asen-Palmer 


Present work 


L (mm) 


2.4 


3.8 


D = 4.4 and L = 20 


F 


0.8 


0.8 


none 


A (s 3 ) 


2.410" 44 


1.778610" 44 


3.510" 45 


B TO (K- 3 ) 


1.010 -11 


1.510 -11 


7.310~ n 


Btu (s) 


5.010" 18 


4.510" 18 


0.8910~ 18 


Bl (K- 3 ) 


6.910" 24 


9.010" 24 


8.610" 24 



TABLE III: Parameters for germanium relaxation times 



p 


0.8 


0.81 


0.82 


0.83 


0.8 


0.86 


0.9 


0.94 


0.97 


1 


T(K) 


10 


20 


40 


60 


80 


100 


125 


150 


175 


>200 



TABLE IV: Ratio of specular to diffuse reflection used in figure [13] for silicon nanowires 
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